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Abstract 



(1) The temperature dependence of the specific heat for a marginal Fermi liquid has been 
calculated . We showed that the expected TlnT correction is characteristic for the low tem- 
perature domain. The high temperature domain has a supplementary correction. The results 
are in agreement with the non-Fermi behavior of some metallic systems in the low temperature 
domain. 

(2) We calculated the self-energy at T = for a two dimensional fermionic system with 
hyperbolic dispersion. The existence of the saddle points in the energy gives rise to a marginal 
behavior, a result which has been obtained by numerical calculations. 

(3) We present a simple demonstration that the two-dimensional fermionic system with the 
energy ek = kxky has a non-Fermi behavior. The calculation of the wave function renormalization 
constant Z were performed using the "poor man's renormalization" method and we have showed 
that Z ^ with the infrared cut-off A as Z{A) = A'' where is a constant. 

(4) The electronic self-energy due to the electron-spin interaction is calculated using the 
one-loop approximation for the two dimensional system and quasi-two dimensional (anisotropic) 
model. We analyzed the relevance of the diffusive modes and the temperature dependence of the 
magnetic correlation length for a possible temperature dependence of the pseudogap. 

(5) We study the influence of the amplitude fluctuations of a non-Fermi superconductor on 
the energy spectrum of the two-dimensional Anderson non-Fermi system. The classical fluctua- 
tions give a temperature dependence in the pseudogap induced in the fermionic excitations. 

(6) Using the field-theoretical methods we studied the evolution from BCS description of a 
non-Fermi superconductor to that of Bose-Einstein condensation (EEC) in one loop approxima- 
tion. We showed that the repulsive interaction between composite bosons is determined by the 
exponent a of the Anderson propagator in a two dimensional model. For a ^ the crossover is 
also continuous and for a = we obtain the case of the Fermi liquid. 

(7) Using the renormalization group approach proposed by Millis for the itinerant electron 
systems we calculated the specific heat coefficient 7(r) for the magnetic fiuctuations with sus- 
ceptibility ~ 5" + + f{q) near a Lifshitz point. The constant value for a = 4/5 and 
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the logarithmic temperature dependence, specific heat for the non-Fermi behavior, have been 
obtained in agreement with the experimental data. 
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Chapter 1 

Non-Fermi Behavior and the Fermi 
Surface 



One of the most important problems which was generated by the discovery of cupratcs supercon- 
ductors is the nature of the normal phase. In fact the main question is if a conventional Landau 
theory of the Fermi liquid (FL) can describe this phase, or as was suggested by the experimental 
date, this phase is in fact a special metal describe by a non- Fermi liquid (NFL). 

The first model of a NFL was presented by Varma et.al.[l] and is well known as the marginal- 
Fermi liquid (MFL) and developed by Kotliar et.al.[2] and Littlewood et.al.[3]. This model is 
phenomenological but due to its importance for the explanation of the experimental data much 
effort has been done to give a microscopic foundation for the occurrence of the MFL behavior. 
The main idea of the MFL model is the occurrence of a linear energy dependence in the imaginary 
part of the self-energy due to the anomalous frequency and temperature dependence of the 
density-density ( spin-spin ) response function. This particular form of the self-energy gives rise 
to the pair-breaking effects, studied by Horbach et.al.[4] and Bendle et.al.[5]. 

This particular behavior of the quasiparticle spectrum is essential for the thermodynamics of 
the system, and the specific heat is one of the most important physical quantities that has to be 
affected as mentioned by Reizer[6] for the electrons interacting with a transverse electromagnetic 
field. 

For a three dimensional (3D) Fermi liquid, the imaginary part of the self-energy has the form 
/mE (lo) ^ and gives a logarithmic correction to the specific heat Cy = + 73£)T^lnT. 
This behavior has been obtained in the framework of the Landau theory, but this correction 
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cannot describe the low-temperature behavior of the ^He which is considered as the standard 
Fermi hquid. 

Anderson [7] showed that the low temperature dependence of the specific heat for ^He can 
be fitted by a T In T dependence and suggested that this dependence is given by coupling of the 
fermionic qausiparticles to the collectives modes. Balian and Fredkin[8] and Berk and Schriffcr[9] 
showed that the coupling between the fermionic quasiparticles and the triplet paramagnon gives 
a T^lnT correction in the specific heat. 

The correction to the linear dependence of the specific heat for a Fermi liquid has been 
calculated by Amit et.al.[10] in terms of Landau parameters. 

Carneiro and Petrick[ll] performed a microscopic calculation of the thermodynamic potential 
in terms of fully renormalized single-particle propagator considering the effect of the quasiparticle 
lifetime. This method can be successfully applied for the calculation of the specific heat of the 
MFL [12] 

1.1 Thermodynamic Properties of the MFL 
1.1.1 Thermodynamic Potential and Entropy 

In this section we calculate the specific heat of a MFL using a similar method with that of 
Carneiro and Petrick [11]. The thermodynamic potential f2 is given by the fully renormalized 
single-particle propagator G (p, u) as: 

O = T ^ e-n^ (p^ ^^)] _ ^ (p^^^) Q (p^ (11) 

n 

where the propagator G (p, Wn) is given by Dyson equation: 

(p, a;„) = Go ' (p, a;„) - S (p, a;„) (1.2) 

E (p, Wn) being the self-energy. 

The thermodynamic potential fl has the property that it is stationary under variation of G 
at fixed Go : 



Jg) 



= (1.3) 

Go 



and shows that we can neglect the temperature dependence of the spectral density A (p,a;). 
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A (p, a;) = -2ImG (p, a;) 

in the expression for the entropy S (T) given by: 



(1.4) 



Using: 



and: 



G (p, u) = ReG (p, u>) - -A (p, u) 



(1.5) 



(1.6) 



E(p,a;) = i?eS (p, cu) - -F (p, u;) 



expression (1.5) becomes 



da; On (w) 



2Im In 



-G-^ (p,w + iO+) 



-Ey_ ^^i?eE(p,.)A(p,a;) 

-Ey_ ^^i?eG(p,a.)r(p,a.) 
where n (u) is the Fermi-Dirac distribution. This equation has been written as: 



(1.7) 



(1.8) 



S{T) = Sdq{T) + S (T) 
where Sdq (T) is the dynamical contribution expressed by: 



(1.9) 



'(r) = E/_ 



+°° doo dn 



CO 



27r dT 



(2Im In [-G-^ (p, a; + iO+)] - ReG (p, w) T (p, a;„)) (1.10) 



The second contribution has been identified as the contribution of terms that has vanishing 
energy denominators.Equation (1.8) can be written as: 



S{T) = So{T) + Si{T) + S2{T) 



(1.11) 



where: 
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So{T) = -Y: [n{el) ln(n(£0 )) + (1 - n(4)) ln(l - n(4))] (1.12) 
p 

and: 

^i(^) =-i: /_r^^^^^(P'-)^(P'-) (1-13) 

For 52 (T) we get: 



/- 



+°° duj dn (w) 



5(p,a;) 

arctan B (p, w) + 



1 + S2(p,a;) 



(1.14) 



where: 



and £p is the solution of the equation : 



i?eG~^ (p,ep) = (1.16) 

If the damping of the quasiparticles is small, Sq{T) is the entropy of the Fermi hquid and gives 
a linear contribution in the specific heat of the form: 

(T) = tT (1.17) 

In fact such a result will be obtained if we consider that in eq.(l.lO) G(p,a;) = Go (p,w) where 
Go(p,a;) is the propagator for a FL, and in this approximation -S'i(r) and S'2(T) appears as 
corrections to the linear dependence. These results will be used for the calculation of the specific 
heat in the following section. 

1.1.2 Specific Heat 

A marginal Fermi liquid can be characterized by the functional form of the frequency and 
temperature dependence of the single particle self-energy [1]: 

^ , \uj\n^ -i^XT uj<T 

EH= ' (1-18) 

Aw In ;f - if Aw T <uj <Uc 
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where A is a coupling constant and ujc is an ultraviolet cut-off.According to Varma and co- 
workers [1-3] the normal phase of high-Tc oxides is an MFL that can be understood on the basis 
of an MFL hypothesis i.e. the total electronic polarizability has a contribution V representing 
an electronic excitation with the property: 

ImV{p,u) = { ^ (1.19) 

-N{fd)sgn{uj) \uj\ > T 

for the leading frequency contribution, where N (0) is the bare electronic density of states. All 
the universal anomalies as well as the appearance of a 'Fermi surface' within the bound states 
can be understood from this single hypothesis. 

From eq. (1.18) a contribution to the resistivity proportional with T is obtained [1]. The 
nuclear relaxation time Tf ^(T), tunneling conductance g{y), optical conductivity u{uj), fit well 
the experimental data. 

The electronic contribution to Cy, in the normal state cannot be extracted from the data 
with high accuracy because the contribution of phonons is very large above the high Tc of these 
materials. k{T) in the normal state is observed to be nearly temperature independent [13]. If 
the self-energy from the eq.(1.18) is used to calculate the entropy, logarithmic correction to the 
linear-dependence are found. The observed n/Ta ratio allowed a weak correction. 

The spectral density A (p,a;) can be calculated from the eq.(1.4) and (1.18) as: 

A (p, a; = ^ ; 1.20 

The contribution given by S\{T) has the form: 

I 2Ar(0)T fFl'i) -F(0)) In^ w < T 

S,{T) = \ \ \' ^ (1.21) 

2N{0)T [F (^) - F (i)) In f u;<T<u, 



where N{0) is the density of states and F{x) is defined as: 

F{x) = 4x-x'^ + x'^tanhx + 2dJCt&n{e'') + 2ReL2{e'') (1.22) 



and: 



L2{x)= / — ln|l-t| (1.23) 

JO t 

The contribution given by S2{T) was calculated using eqs. (1-14), (1.18) as: 
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S2{T) = { 



2Ar(0)T(F(i) -F(0)) (1 + Aln^^) u<T 
2iV(0)r(F(f^) (1 + Alnf ) uj<T<uj, 
Using the results for the entropy we can calculate the specific heat as: 



(1.24) 



Cv{T) = T^^ (1.25) 

In the following we will consider only the contribution given by the MFL, So{T) giving a linear 
contribution specific to the FL. 

^Prom eqs. (1.21) and (1.24) we get for < a; < T : 

Cy = 2iV(0)r (^F -F{0)j (^l + 2A(^lni^-l)) (1.26) 
For the limit T < cj < Uc we get: 



Cv = 2JV(0)r(F(g,)-F(i))(l+2A(ln!^-l))- (1.27) 



(g)r(g)(.(^) -.(!))( 



1 + 2A ( In ^ - 1 



where F'{x) = j cosh^(x). 

From eq.(1.26) we can see that for the MFL the temperature dependence has a correction 
of the form TlnT but in the dominant region T < a; < o^c a more complicated temperature 
dependence is added to this correction. 



1.2 Marginal Behavior of a System with Saddle Points in the 
Energy 

The non-Fermi behavior of electrons with hyperbolic dispersion = ^x^y has been predicted 
by News et.al. [14] in connection with an electronic mechanism of superconductivity in the 
high critical temperature superconductors. For such a two-dimensional (2D) system the density 
of states shows a logarithmic energy dependence, and for a Fermi energy Ep near the saddle 
point (Fc) the critical temperature for the superconducting state has been calculated, and is 
very sensitive to the new energy scale E* given by the density-density excitation. The results 
from [14] are not realistic because only the electron-hole channel has been considered [15], 
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but the results remain important for the explanation of the anomalous properties of cuprate 
superconductors in the normal state. In order to explain the physical properties of the normal 
state of these superconductors we have to mention that the difficulties are due to the study of a 
2D interacting Fermi system and one of the simplest model is to treat it as a Luttinger liquid, 
which exist for ID Fermi liquid. However for the (quasi)2D systems which are more realistic for 
cuprates superconductors, there is no general demonstration for a MFL behavior. As it will be 
showed the form of the Fermi surface is essential for the marginal behavior. A more realistic 
model seems to be the model in which the Landau concept for the quasiparticles is still valid 
and the deviations are due to the anomalous scattering of electrons on the overdamped low 
energy excitations. Millis, Monien and Pines [15] developed such a phenomenological form of 
the dynamic susceptibility x(q, a;) which describe the normal and the superconducting phase, 
the occurrence of the magnetic pseudogap and the influence of the impurities. The pairing has 
also considered in this model, but the Migdal Theorem discussed by Grosu and Crisan [16] 
is not valid in the "hot spots" of the Fermi surface how was recently pointed by Amin and 
Stamp [17]. Another interesting model has been proposed by Bernard et.al. [18] and is based on 
the interaction between the electrons with a critical bosonic mode. This interaction lead to a 
non-Fermi behavior (even to a Varma et.al. [1] marginal behavior) and Crisan and Tataru [19] 
showed that this is possible even in 3D case. The importance of the van-Hove singularities in 
the density of states in the 2D Fermi systems for the cuprate superconductors has been pointed 
by Priedel [20], Labbe and Bok [21] and in a systematic study by the IBM group [22]. In one of 
their papers [23] it was suggested that the saddle points singularities (called "hot spots") of the 
Fermi surface are responsible for the deviations from the usual metallic state. In the following we 
calculate the self energy of 2D model with dispersion = kxky/2m at T = 0. Such a calculation 
was also performed in [24] but the integrals were numerically evaluated. We will give an analytic 
calculation which show the "marginal behavior" of 2D systems. 

1.2.1 Self-Energy of the 2D Fermi Systems 

We consider the self-energy of a 2D electronic system interacting with the electron-hole excita- 
tions which has the polarizability given in [23-24]: 

/mx(q,z) = — [|z + £k|-k-£k|] (1-28) 
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Rex (q, z) = 



1 

2^ 



In 








+ 2 








Z - Ek 





where £k = kxky/2m . The self-energy at T = is given by the general equations: 

ImT, (p, a;) = / - — / dzlmx (p — k, z) /mG (k, z + a;) 

TT J (27r)^ Jo 



-ReS (p, w) = / / dzRex (p - k, z) /mG (k, z + 

TT J (27rj^ Jo 



(1.29) 



(1.30) 
(1.31) 



where u is the electron-electron interaction which generate the bosonic excitation. We mention 
that the dispersion 



(1.32) 



called the hyperbolic dispersion in [24] is equivalent to the simple case £k = k^ky/im. In order 
to perform the calculations using eqs.(1.28) and (1.29) we will perform the transformation: 



and the dispersion has the form 



kx = ^[\ + ^ 



I. - ^ 

fey - 2 



A- 



A 



2m 



(1.33) 
(1.34) 



(1.35) 



1.2.2 Imaginary Pcirt of the Self-Energy 

Equation (1.30) can be transformed using: 



ImG(k, CO + z) = —7r6{z + to — Ek) 



as: 



(1.36) 



7mS (p,a;) = ^ / kdk [| 

TT'^ Jo J-oo A £p_k 



£k — o; -I- £p-k| — kk — <^ — £p-k|. 



(1.37) 
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Let us calculate this integral first in the limit lj <^ e^- In this case eq.(1.37) becomes: 



y2 POO i-c 

ImTi (p, w) = — ^ / kdk [e-^ - co] 

TT^ Jo J-i 



2m 



A {pi +kl- 2p^k^) - + kl- 2pyky) 

and if we take the p = (p, 0) this equation becomes: 



2mu 

ImTi (p, w) = — 5 — / kdk [sk — ^] 

TT^P Jo J-c 

and performing the integration over A we get: 



dX 



X2 _ (p2 _ jt2)/(pA;) + 1 



(1.38) 



(1.39) 



ImT, (p, a;) 



2mu 

TT^P Jo 

2pk 



kdk 



2m 



■ to 



, 2pk — p^ — k"^ 
arcth 7, — h 



2pk n"^ + k'^ 

rarcth- 



p2 _ ^2 p2 _ f^2 p2 _ p'^ — k'^ 

In the approximation arcth{x) = x we calculate the first contribution in eq.(1.40) as: 



/mSi(p,a;) = ^ — 

■TV" n 



rOD 


"k2 


/ dk 


u; 


Jo 


2m 



2pk 



(p + k)2 



— LU 



(1.40) 



(1.41) 



and for <C p eq.(1.41) will be written as: 

2mu^ roo 

/mSi(p,a;) = 3-^/ kdk 

TT^p-* io 2m 

where we neglected the contribution given by fc^ + 2pk <C p^ in the denominator of eq.(1.41). 
From eq.(1.42) we obtain: 



(1.42) 



/mSi(p,a;) 



mu 



2m 



(1.43) 



In the same approximation we calculate the contribution from the second term of eq. (1-40) 



and: 



/mE2(p, w) 



2mu 



2m 



— CO 



(1.44) 



Prom these equations the imaginary part of the self-energy becomes: 

3mu^ 



/mS(p, uj) 



TV" 



P 

2m 



(1.45) 
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and for p = pp we get: 



/mS(p, u) = -3N(0) (w - Ef) 



(1.46) 



where A''(0) = m/7r and Ep = p\/2m 
1.2.3 Real Part of the Self-Energy 

The general expression for i2eS(p, uj) in the second order of perturbation theory has the general 
form: 



ReT.{p,Lv) = - — J dzRex(p - k,z)ImG(k,z + 



(1.47) 



where ImG will be considered as (5-function and i?ex(k, z) is given by eq.(1.29). If we perform 
the integral over z eq.(1.47) becomes: 

i2 roo 



In 



4Ecep_k 2 



2 '^p-k 



£k 

£p-k 



In 



Ek — w + ep_k 



+ 2 



(ek - uj) 
6{w - £k +'^) 



(1.48) 



£k - - £p-k 

where is a frequency cut-off. The integral containing the first two terms will be transformed 
as: 

,2 .oo , 



i?eS(p,a;) 



u 



, , kdk / — 



1-^^i^ln 
Sp-k 



Ek - W 



Sp-k 



-1 + 



"''-'"in 


^^.k-c. 


In'P-'^l 


Ep-k 


£p-k 


4^c 



Using the expansion: 



(1.49) 



ln(l + x) = X — 



we get: 



iieE,(p,u,) = -^/^°°kdk/ 



dA 
T 



gk 

£p-k 



+ In 



4E. 



£p-k 



In order to perform the integration over A we calculate 



(1.50) 



(1.51) 
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£p_k =p^ + k^- 2pk 



(1.52) 



where we take (2m = 1) and (1.51) becomes: 

■■2 foo 



i2eSi(p,a;) = — 



u 



kdk 



47r3Jo |(p-k)(p + k] 



r0(tU + o; — £k) 



(1.53) 



which gives: 



_ p2 J 



Re'Ei{p,oj) 



u 

8^ 



(w + uj) + (ep - w) In 



w + OJ 



- 1 



u 

2^ 



In 



+ a; 



- 1 



The last term of eq. (1.49) is: 



i?eE2(p, w) = -7-^ / kdk / —6{zu - eu + io 
^vr^ Jo J A 



we have to consider k"^ < w + lo and using: 



(1.54) 



:i.55) 



(1.56) 



A 



(1.57) 



k y/zU + LU 

we obtain that i?eS2(p,w) vanished for p — pp — > 0. The final form of the real part of the 
self-energy is: 



ReT,{p,Lo) = 



u 



2 r 



87r3 



{w + + (Ep — cu) In 



Ep 



u 

2^ 



In 



Ep 



(1.58) 



We can see that the divergence in the real part has the form: 



ReT,{p,uj) = (a; - Ep)ln 



Ep 



(1.59) 



These results have been obtained in the one loop approximation. In the following we will 
show that this model present a non-Fermi behavior using the renormalization group method. 
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1.3 "Poor Man's Renormalization" for a System with Saddle 
Points in the Energy 



1.3.1 Model 

Using the "poor man's renormalization " we will present a direct demonstration that the saddle 
points can give the non-Fermi liquid behavior. 

We consider, following Newns et.al. [14] a 2D fermionic system with the energy: 

£k = kxky (1.60) 

where \kx\ < kc;\ky\ < kc, kc being the cut-off. If the Fermi level is closed to the van-Hove 
singularity the nearly 2D fermionic system present two important differences to the normal 
metal. The first one is that the phase space for scattering is much less restricted that in the 
conventional metal and the lifetime of the quasiparticles becomes of the order of their energy. 
Such a fermionic system becomes a MFL. The second point is that the presence of the van- 
Hove singularity close to the Fermi energy is destroying the nesting of the Fermi surface which 
gives magnetic instabilities which are competing the Cooper instability. The effective interaction 
between the quasiparticles has been calculated in ref.[14] as: 



Veff = r (1.61) 

n(q,a;) 



where the polarizability n(q, w) is given by: 



n(q,z) = ^ 



(1.62) 



where / is a fraction of the area of the Brillouin zone and D = k"^. The function P{z,x) is 



a iicii^tujii ui uiic ciicci ui uiic xjiiiiuuiii ziuiic ciiiu ±j — h'^ 

defined as: 



P{z,x)=g{l)-g{-l)-g[-\ + g (--\ - g U --\ +g(-l+'L\ (1.63) 



xJ \ xJ \ xJ \ X 



where 



g(u) = {z + u)\tl\z + u\ (1.64) 
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Using in eq. (1.63) the expansion: 



Inx ^ 2 



we can approximate (1.62) by the expression: 



x-1 
'x + 1 



P{z,x) 



4z 4 
+ - 



x'^ X x^ 



the effective potential has the form 



(1.66) 



(1.67) 



f 

Vr 



;i.68) 



where Vq is given asVo = Ze and A are constants, and Vc is the saturation value of Veff{q,, v) 
in the high frequency region and E* = [fVc — ZeVof /D . 



1.3.2 Self-Energy and the Renormalization Constant Z 

We start with the general expression for the self-energy: 



E(k,a;) = iy J ^Veff(q,i^)Go(k-q,a;-i/) (1.69) 



d^q f du^ 

where Go is the Green function for the free fermionic particles and Vejf is given by eq.(1.61) 
The integral over will be written as: 

d// v. 



f 

/(k-q,a;) = i / 

J — 1 



27r LU — U — 6k-q + iS 



+ if^ (1.70) 



"0 dv I 1 

2tT U — V — £k-q + i5 



where uq = ^jE*e^ and is the matching energy. In order to calculate Z the renormalization 
will be performed only on the real part of the self-energy S^(k,u;) = ReS(k,u;) and we will be 
interested only in the contribution: 
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Rel (k - q, a;) = -Vc9 {yo -\oj - £k-q|) 

+ (^1 - K) e (yo -\^- £k-q|) (1.71) 

+A l^~!!l-ql g(z.o-|a;-£k-q|) 

where ^(x) is the step function. 

Using this result the real part of the self-energy (1.69) is: 

S>,a;) = Ef(k,a;)+A j ^^^^^^J^e {uo - \u - e^^^\) (1.72) 

where S{^(k, cj) gives no contribution to the renormalizcd part. In order to calculate the renor- 
malization constant Z we will use the "poor man's renormalization " taking Sq as shown in 
Fig.1.1. 

If we introduce near kc the cut-off A from eq. (1.72) we calculate dS^ as: 

dE«(A,k,a;) = 4r ^ ^h^£^Szsl .f,^ _ MU+ 
2 JK-dK 27r Jo 27r y " Qx J Qx 

2 JK-dK 27r ^0 27r y Qy J Qy 

which gives 

d^^{A,k,u}) = -A—uj + AvFk— (1.74) 

where: 



dE^(A, k, a;) = 5]^(A, k, w) - S^(A - dA, k, w) (1.75) 

In order to obtain the general equation of the renormalized quantities we write the propagator 
G(A, k,a;) as: 

G(A, k, a;) = ^^^^j (1.76) 

^ ^ a;-VF(A)k-hi7 

and: 
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Figure 1.1: The domain of integration for the calculation of self-energy 



G(A,k,a;) ^ _ ^^^^^^ + 5]R(A,k,a;) + i5]i(A,k,a;) ^^''^^^ 
where we approximated e{k) = vpk , k being the momentum measured from the Fermi surface 
momentum kp- From these two equations we get: 



dZ{A) - VFiA)k) - dZ (A) S^(A, k, - Z(A)55]^(A, k, w) = -dvF(A)k (1.78) 

and in this equation the second term can be neglected since it given no contribution to the 

renormalized quantities and it is in fact a constant which can be absorbed in u. Therefore 
eq.(1.78) becomes: 

{oj - vpk) - dT.^{K) = -dvF{A)k (1.79) 

Z(A) 

on the other hand if dJ2^W has the general form: 
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- dS^(A) = Ci^u + C2VFk^ (1.80) 
A A 



we get the general equations: 

1 dZ{A) _ Ci 
Z{A) dA ~ T 

and 

1 dvFiA) Ci - C2 



(1.81) 



vpiA) dA 
iProm the eqs.(1.81)-(1.82) we get: 



(1.82) 



vf{A) = (1.83) 

and 

Z{A) = A^ (1.84) 

Equation (1.83) shows that the Fermi velocity is constant and Eq.(1.84) shows the important 
result : 



with an exponent 



fk 

and using the relation D = we obtain the exponent 



hm^Z(A) = (1.85) 



^ ZeA^/D 

C= „, (1.86) 



C = ^ (1.87) 



/ 

which depend only on the constants contained in the effective potential 



1.4 Summary of the Results 

•In the first part of this chapter we performed an analytic calculation of the electronic specific 
heat of MFL and we showed that the T In ^ correction predicted in [3] appears only in the 
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domain < a; < T. For energies in the domain T < lo < lOc & supplementary correction 
depending also on In ^ appears. 

Following the same method in the nonmarginal case the correction in Cy is In T if the 
electron-electron or electron-phonon interactions are considered for 3D FL model [26]. This 
behavior makes very difficult the separation of MFL contribution as mentioned by Varma et.al.[l- 
2]. Recently Coffey and Bedell [27] showed that a 2D Fermi liquid has analytic correction in T of 
the form dCy^ = ^2dT'^ and there is no breakdown of the Fermi liquid in 2D. On the other hand 
we mention that for the uranium compound i7o.2^o.8-Pc?3 [28] was found a TlnT correction in 
the specific heat, a result that shows the existence of the non-Fermi behavior in the real systems. 

•In the second part of this chapter we showed analytically that the 2D electronic system with 
hyperbolic dispersion exhibits the marginal behavior ( ImJ^ipF,^^) ~ ^) obtained numerically 
in [24] . The marginal behavior is given by the saddle points of the Fermi surface as was predicted 
in [23]. We have to mention that the marginal phenomenological model is not restricted to the 
2D Fermi systems but the cuprates superconductors are (quasi)2D systems. 

• Using the "poor man's renormalization " method we showed in the last part of this chapter 
the occurrence of the MFL behavior for 2D fermionic systems with saddle points in the energy. 
Finally we mention that the NFL behavior is typically for the ID fermionic systems (Luttinger 
liquid ). The 2D fermionic systems with saddle points in the energy is very realistic for the 
cuprates superconductors, but cannot be treated analytically without renormalization group 
methods. More than that even the NFL models based on the coupling to a gauge field [29] or 
on the superlong range interaction [30] has been treated by this method. 
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Chapter 2 



Anderson non-Fermi Model 

The microscopic description of the superconducting state in cuprates superconductors is a very 
difficult problem because at the present time it is generally accepted that in the normal state the 
elementary excitations are not describe by the Fermi liquid theory. However using the BCS-like 
pairing model the Gorgov equations have been applied to describe the superconducting state in 
the hypothesis that the normal state is a non- Fermi liquid described by the Anderson model [1]. 
The superconducting state properties have been discussed by different authors [2-8] and even 
if these descriptions are phenomenological, it can be a valid starting point for a microscopic 
model. Recent experimental data (ARPES) showed that these materials presents even more 
remarkable deviations from the Fermi-liquid behavior due to the occurrence of the pseudogap 
at the Fermi-surface. 

The occurrence of the pseudogap has been explained using different concepts such as the 
spin fluctuation [9], preformed pairs [10], S0(5) symmetry [11], spin-charge separation [12], and 
the fluctuation of the order parameter induced pseudogap [13]. 

2.1 The Model 

The non-Fermi liquid behavior of the normal state for the cuprates superconductors proposed 
by Anderson [1] was developed by different authors in order to describe the superconducting 
state in the framework of the BCS theory. In the normal state the electrons are describe by the 
Green function: 
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G (k, lUa) 




(2.1) 



1-a 



where uJc is a cutoff energy and < a < 1 and g{c() = 7ra/2sin(^) and cp = — 7ra/2. The 
Green function given by eq.(2.1) was given by Abrahams [14] and used by Yin and Chakravarty 
[15] to study the non-Fermi hquid superconductors. However, the time reverse symmetry is 
violated unless <p = —Tra/2. This is because of the invariance with respect to a; — —u and 
ffk — £^k- When = 0, the spectral function is similar to that of Luttinger liquid without 
spin-charge separation in d = 1 [16] but is very different when 7^ 0. To preserve the equal- 
time anticommutation relation of the fermions, the spectral function must satisfy the sum rule 
dLjA{\<., Lo) = 1, from which relation it can be found the value of g{a) given above. With the 
above value of g{a) the anticommutation relations will be satisfied if |ek| < ^c, which is precisely 
the regime in which a scaling theory is appropriate. The density of states can be calculated from 
the spectral function [15]. In contrast to the Luttinger liquid the density of states does not 
vanished as lo ^ 0. In fact it is unchanged from the Fermi liquid value as a 0. It would agree 
exactly with the Fermi liquid value, if uJc is chosen to be VF(the band width) 
The free action for the model is: 



where the Green function G is given by eq.(2.1). If we follow Shankar [17] we see that the four 
fermions interaction is irrelevant. The spectral anomaly is more stable than the Fermi liquid. In 
fact in the weak coupling regime it does not even allow a superconducting instability [15]. The 
coupling has to reach a treashold before superconducting instability occurs. 

2.2 A Novel Feature of the non- Fermi Behavior. The Occurrence 
of the Pseudogap. 

The recent angle-resolved photoemission spectroscopy (ARPES) experiments confirmed the oc- 
currence of pseudogap in the density of states of the electronic excitations from the normal and 
superconducting phase of the cuprate superconductors. 

Many theoretical approaches and Monte-Carlo simulations [18-19] of strongly correlated elec- 
tron systems has been performed in order to explain the experimental data. Some common fea- 




(2.2) 
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ture are observed in many of these models, the most important being the fact that the origin 
of the pseudogap is the interaction between the electronic excitations and fluctuations. This 
idea leads to the hypothesis that the pseudogap appears in the proximity of a quantum phase 
transition (QPT) [20-21] 

A consistent picture is emerging from the study of the Hubbard model in 2D and 3D using 
Monte-Carlo calculations and the analytically non-perturbativc many-body approach perform 
by Tremblay group. The main results obtain by authors in Ref. [22-25] is the occurrence of a 
temperature independent pseudogap for a 2D Hubbard model taking the interaction between 
the electrons and the renormalized classical fluctuations. In this regime the coherence length is 
increasing with decreasing of the temperature and a good approach for the model is a QPT. 

On the other hand the spin fluctuations are present in both normal and superconducting 
phase of the cuprates. The Fermi surface evolution and the temperature dependence of the 
pseudogap have been studied by Chubukov [26-27] using the phenomenological model proposed 
by Millis Monien and Pines [28] (MMP). 

In this section we will perform an analytic calculation of the pseudogap given by the inter- 
action between electrons and magnetic fluctuations. Starting such a problem we expect that a 
tractable approach will be only one-loop diagram for the self-energy. This approximation does 
not take the vertex corrections which are important because in such a system probably the 
Migdal theorem is not valid. However, even such a simple approximation can gives us an idea if 
the interaction between electrons and spin fluctuations are good candidate to explain pseudogap. 

The self-energy is essential determined in the model by the dynamic susceptibility i^) 
which will be discussed below. Such a discussion is also usefully because at the present time 
there are models with x(k, w) pure diffusive [29] or x(k, a;) describing overdamped fluctuations. 
Recently it was presented a model with both contributions in x(k, w) [30]. 

2.2.1 Dynamic Susceptibility 

We will consider that the electronic system interacts with magnetic excitations near a quantum 
critical point (QCP) determined by the condition J(Q) = — Xo ""^(Q)' where J(Q) is the magnetic 
interaction and Xo ^(Q) is the bare susceptibility in Q = (=b7r, zLtt). The general form for x(q, cu) 
describing the overdamped antiferromagnetic fluctuations is: 
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where a is of the order of the lattice spacing, 5 is a parameter characterizing the distance 
from the QCP and 7 = -^j^, '^SF being the characteristic energy of the spin fluctuations. The 
renormalization group method (RNG) has been applied to the analysis of the phase transition 
of the itinerant electronic system by Millis [31] taking into consideration the dynamic critical 
exponent z. The correlation time of the order parameter r and the correlation length ^ are 
connected by: 



T = e (2.4) 
which shows that r diverges more than ^. Prom eqs.(2.3) and (2.4) we get: 

e = j (2.5) 

so the 2 = 2, one of the most important predictions obtained in the weak coupling approach. 
At T = (i.e. for the quantum phase transition) the static and dynamics are mixed in contrast 
with the case T > where the dynamic behavior affects the static behavior. 

The presence of the particle-hole continuum leads to an overdamping of the modes associ- 
ated with the order parameter. For the superconductors the critical behavior is in the z = 2 
universality class. However Pines [32] showed that the phase diagram of the cuprates supercon- 
ductors is more complicated. The normal phase was divided in two regions: one with z = 2 
{ujgp ~ ^'~^) which has no pseudogap and the region with z = 1 which present a weak pseu- 
dogap regime ( ojsf ^ = a + hT ) and another presenting a strong pseudogap regime 

( oosF ~ = const ) separated by a line T*{x) where x is the doping. The crossover 

between z = 1 and z = 2 regime was studied by Sachdev et.al. [33] using the scaling analysis in 
the framework of a model. The main result is that fermions do not neccesary overdamped and 
destroy the spin waves. At low damping the crossover is the same as in pure a mode. The differ- 
ence is for a finite damping of fermions in the quantum disorder regime where the energy scale 
can drive a crossover from high temperature quantum critical regime 2; = 1 to low temperature 
quantum critical regime z = 2. 
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In this section we will consider for x(q, a;) the expression: 



I \ Xo(Q) ic. rr-\ 

x(q, = ,^22' ^ (^-^^ 



where A is the gap in the spin wave spectrum. 



2.2.2 Self-Energy of the Electronic Excitations 

In the one loop approximation the self-energy of the 2D electronic excitations interacting with 
a bosonic mode has the form: 

2 f f°° dcj' , coth 1^ - tanh ^ ^ 



(27r)2 27r ^ ' ^ uj + u' - £p+q+Q + iS 

2 

where g is the coupling constant, x(q, a;') the dynamic susceptibility, £p = ^ ~ A* and /x the 
chemical potential. 

There are two limiting cases for the calculation of the self-energy. The Fermi-liquid regime 
appears for ujsf ^ T and the non-Fermi liquid regime in the opposite ( renormalized classi- 
cal ) regime losf ^ T [25]. Perhaps the best known characteristic of a Fermi liquid is that 
Im^(kF,LO, T = 0) r^LO^ and /mS(kF,t^ = 0, T) ~ T^. 

The key to the understanding the Fermi liquid versus non-Fermi liquid regime is in the 
relative width in frequency of x"(q, a;')/a;' versus the width of the combined Bose and Fermi 
functions. In general the function n{u!') + /{u + oj') (n(a;) the Bose-Einstein function, /(w) the 
Fermi -Dirac function) depends on u' on a scale Max{u, T) while far from the phase transition 
the explicit dependence of x"(q, oj')Iuj' is on the scale u>sf ~ S> T. Hence in this case we can 
assume that x"(q, a;')/a;' is constant in the frequency range over which n{oo') + f{uj + oo') differs 
from zero. Hence we can approximate our expression for the imaginary part of the self-energy 
with: 



/I 
— 
TT 



coth h tanh 



2T 2T 



= -5'^(pf) y + (7rT)2J (2.9) 

where: 

A(PP)=/|2M ,2.10) 
J 27r u)' 
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As we can see in this regime /mS(a;) ^ uP, so the system behaves as a Fermi hquid. Near 
an antiferromagnetic phase transition, the spin fluctuation energy is much smaller than the 
temperature. This is the renormalized classical regime. The condition ujsf ^ T means that 
x"(q, 0^') is peaked over a frequency interval a;' <C T much narrow than the interval a;' ~ T over 
which n(w') + f{ijO + Lo') changes. This situation is the opposite of that encountered in the Fermi 
Hquid regime. To evaluate /mS the Fermi factor can now be neglected compared with classical 
limit of the Bose factor ^. Evaluating the integral we obtain the exact value obtained using 
classical fluctuations [25]. 

Neglecting the fermionic contribution and using the expansion: 



ep+q+Q ^ ep+Q + vq cos 9 (2.11) 
and Schwinger representation for the propagator: 

- / dte<'^-'+"^^^ = (2.12) 

i Jo u> — £ + ia 

we can write the general form for the self-energy. Let us introduce the structure factor: 

Sid, T) = ^Imx (q, a;') coth ^ (2.13) 

which describe the behavior of the bosonic fluctuations on different energy scales. Using (2.12) 
and (2.13) the real and imaginary part Z^' = -ReX^ respectively Y," = IfnY^ can be calculated 
from eq. (2.8) as: 

S'(p,a;) = |^y^ dtsin[(u;-ep+Q)t]y^ dqqS(q,T)y^ ^e""*'"' (^-l^) 

S"(p,a;) = -|^y^ dtcos[(u;-£p+Q)t]y^ dqqS(q,T)y^ (^-l^) 
where v = The integral over 6 can be expressed by the Bessel function: 

Ja{x) = ^ dee-''"''^^ (2.16) 
27r Jo 

and eqs.(2.14) and (2.15) becomes: 

ct2 /"OO /■oo 

S'(p, c^) = f- / dt sin [(u; - ep+q) t] / dqqS(q, T) Jo(qvt) (2.17) 

^TT JQ Jo 
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0.2 pco roc 

S"(p,a;) = -|- / dtcos[(a;-£p+Q)t] / dqqS(q, T)Jo(qvt) (2.18) 
^7r JO JO 

These general expressions will be used to calculate the pseudogap in the electronic energy. 
2.2.3 Pseudogap in the Electronic Energy Spectrum 

In this section we calculate the correction to the self-energy due to the interaction of the electrons 
with the spin fluctuations for a; ^ T using eq.(2.7) for dynamic susceptibility. 
^Prom this equation we get: 



Imx{ci,Lj) 



Xo(Q)7^^ 



(,5 + aV-^) +7' 
and the structure factor will be approximated as: 



(2.19) 



5(q,T) 



T7xo(Q) r 



do; 



(2.20) 



where A = 8+a?q^ as we approximated coth ^ ~ The integral in eq. (2.20) can be performed 



analytically in the approximation: 



> 



^SF 



2(<5 + aV) 



as: 



(2.21) 



S(q,T) 



Txo(Q) 



5 + a2q2 

^Prom equations (2.17) and (2.18) using the result: 



(2.22) 



L 



Jo{qvt)qdq 1 



S + a?q^ 



Ko{tv/0 



and: 



(2.23) 



00 2 

dxKoiPx) sin(a2;) = , „ In 
V « + P 



(2.24) 



TT 



poo 

/ dxKo{f3x) cos{ax) = 

Jo \/a^ + 



(2.25) 
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we obtain the real and imaginary part of the self-energy as: 



^(..-ep+gf + (X)2 



a; - Sp+Q + - ep+q) + 



a;-£p+Q - ^(w-ep+q) + 



(2.26) 



S"(P,..) = , ^ (2.27) 

These results are identical with the expressions given by Vilk and Tremblay [25] using for 
the dynamic susceptibility a simple form x(q, 0) = (^~^ + q^)^"*^. The identity of the results is 
given by the fact that the details of dynamics are unimportant as long as it does not change the 
wave vector dependence. 

For the case \u! — £p+q| > | the eqs.(2.26) and (2.27) can be approximated as: 

E'(p,a;) = M!lil!iC^L^£p±^ (2.28) 
6 (jj - Ep+q 

S"(p,u;) = (2.29) 

At this point we mention that the Green function G~^{p,uj) = Go"'^(p,a;) — I](p,a;) will 
present a correct behavior ( G{u!) ~ l/u ) if S(a;) ~ l/u if we have: 

rin^ = c (2.30) 

V 

where c is a constant and Ux is the characteristic energy which is in agreement with the approx- 
imations satisfy Ux ^ T. Prom eq.(2.30) we get: 

and ujx = <C T, define the proximity to the phase transition. This simple procedure can be 
regarded as phenomenology of the Tremblay renormalized classical regime of the fluctuations 
obtained in the self-consistent way, and gives as a main result a temperature independent for 
the pseudogap. 

However, if we consider the phenomenology proposed by Millis, Monien and Pines and take 
z = 1 the pseudogap obtained from eq. (2.28) as: 

Alg = Tln^ (2.32) 
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will present a temperature dependence for the pseudogap. 



2.2.4 Effect of Anisotropy 

The experimental data showed a highly anisotropic gap. The effect of the interplane anisotropy 
has been considered by Preosti et.al. [34]. Following the method developed in previous section 
we will show that there is a temperature dependence of the pseudogap if the anisotropy is 
considered. In this case the structure factor will be taken as: 

Sa {q, qzT) = , , (2.33) 

where = + 7 = ^zKxy and ^^y = iFrom eq. (2.17) we obtain: 



S' (p, u) = ^ r dt sin [(a; - Sp+g) t\ dqq ^ {q, q,, T) Jq {qvt) (2.34) 

27r Jo io io qnz 



I dtsm.\{u - £^+Q)t] I dqq I 

qhz 

In this expression we perform first the integral over variable t, which is different from zero 
only if: 

\lu - ep+ql > vq (2.35) 

and we get the result: 

(P'^) = 2^ ' / ~ ; „ — (2-36) 



•^0 'IBZ ^^^-e^+ctf-vW 



which will be written using approximation (2.35) as: 



(T -i) 

S^(p,a;)= . ^^'^ (2.37) 

(W - £p+Q) 



where: 



A^,(T,7)ro q. 



^LiO) T q%z 



A(C,7) (2.38) 



2^e 

xo(Q) 

and: 



^0 = :r^AvM (2-39) 
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Figure 2.1: Temperature dependence of the pseudogap A(T) for 7 = 0.9 and ^ ^(T) 

exp[iooro/r] 



1 + i&cf 



arctan 



7 



arctan 7^9^ _^ 1 1 + (7^9, 



(2.40) 



We can take a soft cutoff taking Xo(Q) Xo(Q)[l — Qc ^(Q^ + 7Qz)] ^''^d the equation for 



the pseudogap has the form: 



A^„(T,7)ro 



If we define the temperature T* (7) by: 



(2.41) 



^(r*(7)) ^ 00 



(2.42) 



we can show that: 
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T/T„ 



Figure 2.2: Temperature dependence of the pseudogap A(T) for ^ ^(T) = 10 ^ + 10 "^T/Tq,^ 
0.9 respectively 7 = 1. 



limr*(7) = 
7— »o 



(2.43) 



We have to mention that there is an important difference between this result and the result 
from Ref.[34]. In this model ^ is the correlation length of the magnetic fluctuations, differently 
from Ref.[34] where ^ is the superconducting correlation length. The condition (2.42) is equiva- 
lent to the proximity of the magnetic quantum phase transition. Even if this model is realistic 
taking into consideration the magnetic fluctuations the matching between the pseudogap and 
superconducting gap was not solved. An important point of the model is the temperature de- 
pendence of the magnetic correlation length. In Fig.2.1 we present the temperature dependence 
of the pseudogap obtained from eq. (27) and ^~^{T) = 10~^ + 10~^r. Using 7 = 0.9 and 

= 10~^ + lO^^T/To we obtained a dependence of the pseudogap given by eq. (2.38) in Fig. 
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2.2. 



2.3 Electron-Fluctuation Interaction in a non-Fermi Liquid Su- 
perconductor 

In this section we consider the model introduced in [1] and describe in the first section of this 
chapter. We consider that the superconducting state appears due to an attractive interaction 
and is describe by the BCS-hkc order parameter which can be calculated from the Gorkov 
equations. The fluctuations of this parameter can interact with the electrons and the fermionic 
spectrum of the elementary excitations changes. Such an effect has been studied in the BCS 
superconductors by different authors [25], [35] and it was showed that this interaction gives a 
contribution to the density of states for T > Tc which explained the behavior of the tunnelling 
measurements. 

For a superconductor described by the Gorkov-like equations with the normal state described 
by eq. (2.1) the propagator of the fluctuations has the expression: 



where V is the attractive interaction between the electrons and 11 (q, iun) is the polarization 
operator defined as: 

n (q, iu;„) = r E / (p, ia;i) G (q - p, ia;„ - icoi) (2.45) 

where G (p, iwi) is the Green function related to electrons, which in terms of Dyson equation 
has the following for: 

(p, icoi) = Go ' (p, iui) - E (p, iu;i) (2.46) 

where the self-energy is given by: 

S (p, iu,) = -Tj2 f (q, iu;„) G (q - p, ic^„ - ia;i) (2.47) 

Equations (2. 44)- (2. 47) has to be solved self-consistently, but this cannot be done analytically. 
However, in the mode coupling approximation it can be done and we can calculate the new 
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energy of the electronic excitations. We mention that in eq. (2.47) the vertex corrections have 
been neglected in order to developed a simple analytical model. 



2.3.1 Mode Coupling Approximation 

In this approximation we consider first that G (p, iwn) ~ Go (p, \u)n) and from eq. (2.45) we 
define the polarization 



no (q, iu;m) = / 72^*^ ^' ^^m) 



(2.48) 



where: 



5(k,q,ia;^) = (-l)i-"r^ 



-2a 



(2.49) 



We performed the analytical calculation of no(q, iwm) given by eq. (2.48) and from the eq. 
(2.44) the propagator for the order parameter fluctuations has been obtained as: 

rp \2a / rp \ 2a 



L>o' (q,ia;„) = iV(0)A(a)^C(a) 



T \ LOr. UJr. 



+ 



vpq 



where the critical temperature has been obtained as [7] : 



(2.50) 



-■2a 



1 



C{a) 



D(a) 



2« 



A{a)N{iS)V 



and the constants from eqs.(2.50) and (2.51) are: 



(2.51) 



22a 

A{a) = g^ia) sin7r(l — a) 

TT 



(2.52) 



C(a) = r^(a)[l-2i-2°]C(2a) 



D{a) 



r(l - 2a)T{a) 
2ar(l - a) 



(2.53) 
(2.54) 
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20F 



1 - 2 



2-2a 



C(2„-l)(-J 



and: 



B(2- 2a,Q; - 1) /w^ 



2(2a - 1) \uj, 



2a-l 



( ^ 






a, — 








OJc J 





(2.55) 



xi^CP-2„)(^) 



(2.56) 



_ S(3-2a,a-2) /a;z)\2"-2 /T\2 
4(2a - 2) \~J \~J 

B(x, y) = r(x)r(y)/r(x+?/) and T{x) is the Euler function and C,{x) is the Riemann function. 

Using a similar form with one introduced by Schmidt the fluctuation propagator will be written 

as: 



where: 



and: 



Dq^ (q, iwn) = A^(0) \h{a)T{a) + m(a)a;„ + ^^(a, T)q^ 



r(«)H,^) -1 



a{a) 



M 



(a ^ 



(1 - a)A{a) 



rn \ la 

6(a) = A{a)C{pL) ( — 



(2.57) 



(2.58) 

(2.59) 
(2.60) 



(2.61) 



In the approximation S ^ ttT the Green function given by eq. (2.46) will be approximated 
as G = Go + GqSGo and 11 will be modified by 511 also linear in S. Following Ref.[ll] we 
calculated 5Ii in the "box approximation" as: 
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SU = 2T'Y: J ^G§(p,ia;„)G2(-p,-ia;n) J ^D(q,a;„ = 0) (2.62) 

where: 

L>o ' (q, i'^n) = V-^ + U (q, iwn) + SU (q, ia;„) (2.63) 

The box vertex describes the interaction between fluctuations and becomes important only 
in the critical regime in the standard superconductors. In order to calculate (5n(q, ioJn) we 
introduce: 



Bo = ^fY: J ^Gg(p,ia;„)Gg(-p,-ia;„) (2.64) 



NiOy^J (2vr)2 

where A^(0) = m/2n is the density of states. If we used for the electronic Green function eq.(2.1) 
we obtain: 

_ i?(l/2,3/2-2a)[23-^"-l]C(3-4a) u;"^" . 
Bo{T) {nTf-^ ^ ^^-^^^ 

If we introduce T(a) = T{a) + 5Yl/N{Q) the fluctuation propagator given by eq. (2.57) will 

be: 



D-^ (q, iw„) = 6(a)r(a) + ia(a)w„ + (a, T) (2.66) 

where: 

^ 2i3o(T) ^ /■ 1 

- "^"^ = lV(Or^i (2.)^f(a) + e^(a,T)g2 ^ ^ ^ 

If we perform this integral taking the upper limit qm = ^/^{ct,T) from eq.(2.64) we get: 

X .X Bo(T)T , l + ?(a) 

-'°'-^<"> = 2.W(„.r) '°?RfT (^■'"*) 

For realistic parameters ( = lOOiT, a;c = 200K ) the deference f{a) — r(a) becomes 
important only near a critical value of a defined by ^(oc) = 0. In the BCS limit ( a = ) this 
parameter is small and this behavior can be associated with the occurrence of the preformed 
pairs in the domain Tc <T <T*, controlled by a. This behavior is in fact due to the occurrence 
of a pseudogap in the electronic excitations. 
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2.3.2 Electronic Self-Energy 



The self energy due to the interaction between electrons and fluctuations is given by eq.(2.47) 
where D{q,iu!n) is given by eq.(2.66). First we calculate the summation over Matsubara fre- 
quencies Un ■ 

5 = T ^ (q, G (q - p, - icc;i) (2.69) 

S = Ty 1 ij g[a)u, 

N{0){bf + iau;n + eq^){iuJi-i^n + e^)^-'' 

transforming this sum in a contour integral which has a pole at n(q) = — (br+,f^q^)/a and a cut 

line from itoi + ek to oo in the upper semiplane. From eq. (2.59) we can see that a(a) = — \a{a)\ 

and il(q) = — (br + ^^q^)/ |a|. Performing this integral we obtain: 

_ u;-" n(^7(q))g(a)e'W2 sin[vr(l - g)] 

iV(0) Hu;; -ek-!^(q)]i-" iV(0) tt ^ ' 

4 '''[6r-|a|(t + ia;0+eV](i-£k)^-- ^ ^ 

where n(a;) is the Bose-Einstein function and f{x) is the Fermi-Dirac function and £k = k^/'im— 

Ep. The integral from the second contribution in eq.(2.72) will be performed using the expansion: 



fit) = J2 (-1)" exp i-Pim + m (2.73) 

m=0 

and the last term become: 



_ - (-1)"^ [£^ + Qiq)]"/^+^ (I3{rn + l)[(Xq) 
\a\ [/3(m + l)]"/2 



X r(a)W^_,/2,,/2-i/2{/?(m + l)[£k + ^^(q)]} (2.74) 
where the Whittaker function W\^^{z) will be approximated as: 

Wx,„{z) ^ e-'I'^z^ (2.75) 
These results gives for eq.(2.72) the expression: 

^ g;-" n(0(q))g(a)e'W2 . . 

iV(0) [-io;; - £k - ^^(q)]!-" iV(0) 
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sin[7r(l - a)] ^ (-1)"* + m + ^(q) 



TT 



m- 



- \a\ [/3(m + l)]" 



X r(a) exp[/3(m + l)ek] (2.76) 
In the limit k = kp the second term denoted as S2 becomes: 

a;~" sin[7r(l — a)] iuji \a\ + 6r + ^^g^ 



^2 = 



A^(0) TT 



a 



|2 



X r(a)(l - 2i-")C(a)5(a)e-"/2 (2.77) 

and if r Tc, a;/ 0, g this term can be neglected. This approximation is in fact equivalent 
with the physical picture proposed by Vilk and Tremblay [25] in which the occurrence of the 

pseudogap is given by the interaction between the electrons and the classical fluctuations. Indeed, 
in this regime the first term of eq.(2.76) can be written as: 

and the electronic self-energy becomes: 

S (p, u; + iO) ^ -^IgG (k, -iuy) (2.79) 
where we considered ek S> J^(q) and: 

will be approximated as: 

2 _ T fiM qdq 

™ ~ 27riV(0) \a\ Jo (&^+eV)7R 
where qm is the wave number cutoff. From eq. (2.81) we calculate the temperature dependence 
of Al,(T) as: 
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Figure 2.3: The temperature dependence for the pseudogap for £,qj^ ~ kp^ ~ 10, Tc = 100 

The temperature dependence of Apg given by eq.(2.81) using the approximations 'r(a) = t(q;) 
and ln(l + Cq]^) ^ ln(Cg|f) is shown in Fig.2.3 for ^qm ~ kp^ ~ 10,47riV(0)^ ~ 1, and 
Tc = 89K. 

2.4 Field-Theoretical Description of Crossover Between BCS 
and EEC in a non-Fermi Superconductor 

The problem of the crossover from BCS superconducting state to a Bose-Einstein condensate 
(BEC) of local pairs [45-46] becomes very important in the context of high temperature su- 
perconductors (HTSC). While at the present time there is no quantitative microscopic theory 
for the occurrence of the superconducting state in the doped antiferromagnetic materials, it is 
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generally accepted that the superconducting state can be described in terms of a pairing picture. 
The short coherence length ~ 10-20 A) increased the interest for the problem [48-58] because 
it showed that the BCS equations of highly overlapping pairs, or the description in terms of 
composite bosons cannot describe the whole regime between weak and strong coupling. The 
mean field method developed by different authors [47,48,50,58] and solved analytically in two 
and three dimension, and the Ginzburg-Landau description [51,52,57] showed that the evolution 
between the two limits is continues, no singularities during this evolution appearing. The zero 
temperature coherence length in the framework of field-theoretical method has been studied 
in Ref. [56]. The problem of the BCS-BEC crossover in arbitrary dimension d using the field- 
theoretical method has been extensively discussed in Refs. [58-61], where the chemical potential, 
the number of condensed pairs and the repulsive interaction between pairs have been calculating 
using the analogy with the field-theoretical description of superfluidity. In this section we apply 
this method to study the crossover problem for a non-Fermi superconductor described by the 
Anderson model [1] . 



2.4.1 Weak Coupling Limit 



The BCS-like model for the non-Fermi system is described by the Lagrangian 



(2.83) 



where the normal state is describe by the Green function (2.1) 
If we introduce the two-component fermionic field 



/ / ^ 

V'T 



the non-interacting part of the Lagrangian (1) is: 





Lo = I I * 





In order to calculate the partition function: 



(2.84) 



(2.85) 



Z = J D^+D^exp {i j L 



(2.86) 
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we will transform the interaction contribution from the Lagrangian (2.83) as: 



exp (^-iXo Jj+iP+iPi^P^^ = J L>A+DAexp (a+V^V^T + V'^V'J A - ^A+a)) 

(2.87) 

where S^ = JdtjS X as in Ref. [58-61] and A = AqV'j.^t ^ bosonic field. The partition function 
defined by eq. (2.86) will be expressed using eq. (2.87) in a bilinear form as 









A 1 




exp 






j 













Z = J D^+D^ J DA+DAexp (j-A+Aj exp i J\/+ ^ 
Performing the integral over the Grassmann fields the partition function becomes 

Z = J DA+DAexp(^iSeff[A+,A] + yA+A^ 
where Seff[A~^, A] is the one loop effective action,which can be written as: 



(2.88) 



(2.89) 



5e//[A+,A] = -zrrln 



f{a)ipo-epy-- 



J 



(2.90) 



-A 

-A+ f{a){po + ep)'-^ 

where /(a) = uj^g~^{a) and the Tr has been used according to the meaning from Ref. [60]. 

In the mean field approximation the integral from eq. (2.89) can be performed using the 
solution given by the saddle point and for T 7^ the critical temperature Tc is given by the 
eq. (2.51) replacing V with Aq . This expression is valid only in the limit < a < 0.5 and 
a positive critical temperature implies for the coupling constant the condition |Ao| > Ac, with 
Ac = {uJc/u}D)"'/iA{a)/D{a)). Wc have to mention that the critical temperature obtained is 
different from one obtained in Ref. [15] and it is easy to show that it gives the exact BCS result 
in the limit a — > 0. If we consider the effective action as: 

\l-a n \ 



5e//[A+,A] = -irrln 



/(a)(po-£p)'-" 

/(a)(Po + £p)'-" 



-iTr In 1 



A 



and the system as space-time independent the partition function can be written as: 



(2.91) 



Z = Zoexp ( ^A"^A 



(2.92) 
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Zq containing the non-interacting contribution, and we get for the renormahzation couphng 
constant A the expression: 



A Ao + /2(a) J (27r)2 J 27r (^2 _ ,2)^"" ^^'^^^ 



1 1 i f (Pp f dpQ 1 

Using the integral: 

f 1 . ^; ^r(^ - 1/2) 1 

/ = i -1)'^ ^ / ^ 2.94 



we calculate A as 



1 _ 1 _^ 1 
A Aq Ai 



(2.95) 



where 



cos(7r(Q-l))5(l/2, 1/2-a) V^^D/ 
The expression given by eq. (2.96) is positive for a < 0.5. The new coupling constant A has 

to be also negative in order to have superconductivity ( A < ) and this condition is satisfied 

if |Ao| < Ai. If we consider also the condition |Ao| > Ac we get the general condition for the 

coupling constant Aq, Ac < |Ao| < Ai, which is satisfied for < a < 0.5. 

We mention that for the weak coupling limit Aq 0~ the BCS limit studied in Ref. [15] is 

reobtained, but we also showed that the critical constant calculated from the critical temperature 

is Xi{a ^ 0) = 0. 

In the limit Aq — >^ — 00 called the strong coupling limit, we expect an important effect of the 
non- Fermi character of the electrons in the coupling constant. 

2.4.2 Strong Coupling Limit 

In this limit we consider A(x) = A + A{x) and consider the action S'e//[A+, A] obtained from 
eq. (2.90) as: 

Seff (a+, a) = -iTr (1 + GK^ (2.97) 

where: 
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/(a)(po-£p)'-" -A 

-A+ /(a)(po + £p)'-" 



A 



A 
A+ 



which can be written as: 



OO -1 

Seff (A+,A)=-zrr^- 

1=1 



GA 



with: 



A 



-A 
-A+ 



We are interested in quadratic terms in A and we will take the approximation 



(2.98) 



(2.99) 



(2.100) 



(2.101) 



Seff (a+, a) = 5f;^ (0) + (q) (2.102) 

which contains the quadratic contributions. The first term in eq.( 2.102) has the form 

1 



^3/ (0) = l^Tr 



l-o 



A 



A^A+A+ + A^^AA + 2 



+ -iTr 



A 



A 



A 



l-a 



A 



(2.103) 



which will be approximated, taking in the dominator A Ri as: 



5fA(0) = ^^^r 



— f A^A+A+ + A+Va + 2 

2 P{a){pl-el) "V 



A 



A 



(2.104) 



the last term giving no contribution to the renormalized coupling constant. Following the same 

(2) 

approximation we calculated S^j/j^ (q) as: 



Seff (q) - 2^^"/^(a)(Po - ep)i-"(po + 9o - e^^^y--^^'' 



2 /2(a)(po-£p)i-"(po + go-£p+q)i-" 



A+A 



(2.105) 
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^Prom eq. (2.104) and (2.105) we have 



f ^ 

J (27r 



(27r)2 (p2 + m£„)3-4a 



A^A+A+ + A+^AA + 2 



A 



A 



(2.106) 



and 



L(^) (q) 



sin(7r(l — a))Z?(a, a) ^_2q /" d^p 



47r/2(a) 



y (27r 



(27r)2 (p2 _|_ _,_ ^2/4) 



2 //lAl-2a 



sin(7r(l — a))i?(a, a) ^;^_2cj /" d^p 



-m 



47r/2(a) J (27r)2(p2+^£^_^Q^ + g2/4)l-2a 

The integrals from eq. (2.106) and (2.107) can be performed using the formula 



(2.107) 



1 r{N - d/2) 1 



and we obtain 

L(2) = - 



167r2/- 



p (p2 + A^)N (47r)<^/2r(Ar) (^2)iV-d/2 



m 22-4« 5(1/2,3/2 -2a) /-^2-r + x+ ir+2xx n 
^/2(a)l-2a ^ (aA+A+ + A AA + 2 



(2.108) 



A 



A 



m sin(7r(a — 1)) B{oi,a) 
'I67r2/2(a) 2^^ (£» + 9o + gV4m)-'" 



AA+ 



m sin(7r(a — 1)) B{a,a) A"'"A 

^ 167r2/2(a) 2^^ (^a - go + g74m)-'" 
Using the approximation: 



(2.109) 



(£„ - go ± gV4m)''' = 4° + 2a£2a-i / ^ 



4m 



and using the notation 



A 

A+ 



we obtain from eq. (2.109) 



(2.110) 



(2.111) 



sin(7r(l - a))B{a, a)el''-'^ J^r+M* 



167r2/2(a 



(2.112) 
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where 



M 



^0-2^-/^0 



nil, = 2m being the boson mass and /xq the chemical potential 



22-^°5(1/2,3/2-2q!) 



1^0 



P{a) (1 - 2a) sin(7r(l - a))B{a, a) 
the velocity cq of the sound mode is: 



A 



-2a-l 
'a 



CO 



1 



22-'i"5(l/2,3/2 -2a) 



JT^ft /^(a) (1 — 2a) sin(7r(l — a))S(a, a)m 
and the repulsive interaction Ao6 between pairs is: 



A 



.2a- 1 



(2.113) 



(2.114) 



(2.115) 



Ao6 



TT^ 24-4"5(l/2,3/2 - 2a) 



(2.116) 



m (1 — 2a)[sin(7r(l — a))B{a, a)]^ 
We mention that limQ,_>o Ao6(a) = 27r/m a result identical to the result obtained in Ref. [60] 
for the 2D case. 



2.5 Summary of the Results 

• In the first part of this chapter wc performed a simple calculation of the pscudogap due to the 
interaction between electrons and magnetic fluctuations. We showed that for the 2D systems, in 
one loop approximation the self-energy of the electronic excitation can be expressed analytically 
by the structure factor ^(q, T). The expression for the real and imaginary part of the self- 
energy given by eqs. (2.17) and (2.18) are very usefully because their behavior as function of uj 
and q can be easily controlled and predicted if the form of 5'(q, T) is known. This calculation 
make transparent the analytical approach first proposed in Ref. [22-25] using the self consistent 
treatment. 

We considered in the susceptibility x(Q)^) the diffusive contribution of the spin waves and 
we showed that if the condition expressed in eq. (2.21) is satisfied the spin-wave does not change 
the occurrence of the pseudogap predicted by the coupling with low energy spin fluctuations. 
The energy scale has been obtained from the condition G{lo) ~ 1/lo ( 5](a;) ~ l/u ) and in 
fact this condition gives for the coherence length of the magnetic fluctuations an expression 
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of the form ^(T) ~ exp[C/T] used in [25] in the renormahzed classical regime, and which is 
appropriate also for a 2D magnetic system [36]. Recently the RNG equations have been applied 
for the proximity of the Lifshitz point [37] . 

• We also tried to consider the effect of the anisotropy in the magnetic fluctuations and 
following a calculation similar to that of Ref. [34] wc caknilated a pseudogap which is temperature 
dependent. This dependence is even more complicated because in different regimes the coherence 
length ^(r) for the magnetic fluctuations is different. Using different T dependence for ^(T) we 
showed that even for ^(T) ~ exp[C/T] there is a temperature dependence in the pseudogap 
which is not certain from experimental point of view [38]. 

We mention that the key problem of such model remain the matching between the pseudogap, 
which in this case has a magnetic origin, and the superconducting gap which should have a d- 
wave pairing origin. In order to have a picture of this point we propose that the following 
qualitative behavior of the phase diagram. The total gap A(T) is defined as: 



A(r) 



A,(T)+Ap,(r) ,r<T, ^^^^^^^ 

Apg(T) , Tc < T < T* 



where Tc is the superconducting critical temperature and T* is the pseudogap appiaring tempera- 
ture. According to this picture the measured gap below the superconducting critical temperature 
has two contributions, one from the superconducting gap and the other one from the fluctuation 
gap. 

• In the second part of this chapter we showed that a temperature dependent pseudogap 
appears in a non-Fermi liquid superconductor due to the interaction between electrons and the 
fluctuations of the order parameter amplitude. The mode-mode coupling, valid in the weak 
coupling approximation can give relevant results,even for the intermediate coupling studied by 
Levin group [39]. The model recently applied by Norman et.al. [40] can be applied for the spin 
fluctuation model proposed by Chubukov [9] in order to study the temperature dependence of the 
pseudogap. In Ref. [40] the filling in of the pseudogap due to the increment of the temperature 
is given by the broadening in the self-energy and is proportional to T — Tc. A similar broadening 
effect proportional to f{a) was obtained in our model and this can be seen very easily from eq. 
(2.77) if in the electronic Green function we take the limit g = 0. 

Recently such a model for the Fermi liquid superconductor has been studied by Kristoffel and 
Ord [41] and their temperature dependence is different than our result. However, we mention 
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that according to their model these authors have to obtain a result similar to the results given 
in Ref. [25] . The difference is given by the method of performing the integral over q which is not 
correct in Ref. [41]. 

Recently Preosti et.al. [34] generalized the method given in Ref. [25] taking into consideration 
the anisotropy in the dynamic susceptibility due to the interplane pairing. From the temperature 
dependence of the pseudogap shown in Fig. 2.3 one can see that there is a narrow domain of 
temperature where this dependence is in fact in agreement with the recent data from Ref. [41]. 
Increasing the temperature will lead the system to the pure classical regime, where the pseudogap 
is constant, result in agreement with Ref. [25]. 

The model analyzed in Ref. [25] and [41] use a constant coupling between electrons and 
overdamped fluctuations describe by the t-matrix or a x(q, cj) containing an imaginary part. 
Recently Tchernyshyov [42] showed that using a better approximation for the t-matrix the decay 
of the Cooper pairs is negligible and a bosonic propagator without damping of excitations gives 
a pseudogap in the normal state. This idea is also interesting for our model. This simple model 
neglected the vertex corrections in the electronic self-energy, which for the non-Fermi liquid 
superconductor are very singular and has to be considered in the transport theory. 

At the present time it is generally accepted that the pseudogap appears only in the un- 
derdoped cuprates. This fact is reflected in our scenario by the dependence of Apg of a which 
appears as a parameter of the model. However it was showed [44] that the correlation length 
depends on a. and we can attend the overdamped regime by variation of a. This demonstrate 
that such a model is appropriate for the description of the pseudogap in a non-Fermi-liquid 
model. 

•In the last part of this chapter using the field-theoretical methods we studied the crossover 
between BCS and BEC in a non-Fermi liquid. The weak coupling case lead to the same results 
as in the mean field like models. In the strong coupling limit we showed that the pairs form a 
Bose gas with a repulsive coupling constant which is controlled by a 
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Chapter 3 

Renormalization Group Approach of 
Itinerant Electron System near the 
Lifshitz Point 

3.1 Introduction 

The occurrence of the non-Fermi behavior in the systems of ferniions coupled to a critical 
fluctuations mode has been suggested in connection with the neutron experiments [1-2] and 
studied in the framework of many body theory [3-4] in the case of two-dimensional (2D) and 
three dimensional (3D) models. Recent experiments on the heavy fermions systems showed 
also a non-Fermi behavior of these materials at low temperatures and it was associated with the 
proximity of quantum critical point (QCP). The most studied example [5] is CcCuq-x-^Ux where 
at the QCP x = 0.1 the resistivity increases linearly with temperature T over a wide range of T 
and the specific heat C{T) is proportional to Tln(To/T). This behavior has been explained [6] 
by the coupling of 3D fermionic excitations to the 2D critical ferromagnetic fluctuations near the 
QCP. The inelastic neutron scattering measurements performed on this materials [7-8] showed 
the following new points in the behavior of this material 

•The inelastic neutron scattering data can be fitted with a susceptibility of the form: = 
C-^[f{q) + {aT-iuj)°'] where a = 4/5 and not 1 as is predicted by the mean field approximation 

•The quadratic stiffness vanishes, fact which shows that we are dealing with a quantum 
Lifshitz point (QLP) 



50 



•The peaks for x = 0.2 and x = 0.3 can be considered as 2D precursor of 3D order 



•The scaling analysis showed [7] that j{T) = C{T)/T has the form 



7(r) - T^D-l/2)a/2-l 



(3.1) 



which for D = 3 and a = 4/5 gives a temperature independent value. This analysis has 
been performed taking as the most important contribution to x the form containing and 
the q-dependence of the form f{q) = Dq^^ + Cq\ where a = z being the critical exponent 
from the dynamical critical phenomena [9]. In this chapter we will show that using the Hertz 
[10] renormalization group method (RNG) extended for T / by Millis [11] we can obtain the 
ln(To/T) term as a quantum correction to the classical results expressed by eq. (3.1). 

In the next section the renormalization group studies of the Gausian fixed point of magnetic 
transitions in the metallic phase will be use. This treatment of the Gaussian fixed point essentially 
lead to equivalent results to the self-consistent renormaUzation approximation [12]. For the case 
of ferromagnetic system with isotropic Fermi surface we have for the free susceptibility the 
following expansion: 



It should be noted that this expansion is possible only for d > 3. In d = 1 and d = 2 
the susceptibility is singular at |k| = 2kF-, even for an isotropic Fermi surface. Recently the 
nonanalyticy of x(k, io^n) in k for d < 3 in the leading order was considered more seriously 
[14] which led to the conclusion that for the ferromagnetic case, the mean field description is 
incorrect at d = 2 and d = 3 in contrast with the result from eq. (3.1). 

In the antifcrromagnetic case, x(Q + k, io-'n) is similarly expended for small k and LOn if the 
nesting condition is not satisfied and the spatial dimension satisfies d > 3. When the nesting 
condition is satisfied, as in the case of Mott insulator state, this type of simple expansion is not 
possible. 

3.2 The Scaling Equations 

The main idea of our model is contained in a modified form of the susceptibility given in Ref. 




(3.2) 



[7] as: 
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(q, ia;n) = C {f{q) + + {aT^ + la;]") (3.3) 

where f{q) is a smooth function of the wave vector, 5 is the control parameter (pressure, doping) 
measuring the distance from QCP and a is the phenomenological exponent. The real part of 
this expression satisfies for a; = the relation: 



x'-^(q,T)-x'-'(o,r) = c(ar)-" 



(3.4) 



which is identical to the relation satisfied by the form proposed in Ref. [7]. The imaginary part 
of the susceptibility expressed by eq. (3.3) has the form 



x"(a;,T) = r-"5(|, 



E d_ 
T 



(3.5) 



where: 



g{y,x) = C- 



sin 



arctan 



f„.a sin(7ra/2) 



+2/™ cos(7ra/2) 



1/2 



(3.6) 



(1 + x° + y" cos(7ra/2))^ + y2a sin(7ra/2) 
and f{q) = 0. We mention that for 6 = and a = 1 we reobtain the mean field results and it 
can be showed numerically that the scaling function g{y, x) has the same form with the scaling 
function g{y) obtained in Ref. [7]. 

Using these considerations we consider that in the low temperature approximation the in- 
teracting Fermi system can be describe by the effective action: 



where: 



5e// ($) = 5f;/$) + 5^//$) 



(4) 



(3.7) 



si%{^) = VTY: 



(2^ 



si%i^) = uVT' E n ^n)s ( E qi ) M E '-n. 



(3.8) 



(3.9) 



with ujn a bosonic frequency, u > is the coupling constant and D is the stiffness of the 
fluctuations. 

In order to calculate specific heat we will use the Gaussian form of the free energy obtained 
from eq. (3.8) as: 
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dz d^mi dq± , z A sin 6* 

F= — - — i- / -i^coth — arctan- ^ ^ ^ (3-10) 

Jo 2ttJo {2TTf Jq 27r 2T AcosO + qj + Dq\ + q^^ ^ ' 

where Q = arctan(z/(5") and A~°' = ((5^" + z^)^/^. Using the transformations 



to' = b'^/"uj g[j = bqii q'^ = b^/'^q_L 
T' = IP'/'^T S' = b'^/'^S D' = bD 
and following the same procedure as in Ref. [11] we obtain the scaling equations: 



(3.11) 



dT(b) 2 , , 

-TTT = -Tib) (3.12) 
a in a 



'^''^^ -(l-l\uib)-u\n + m2 (3.13) 



d In 6 V 2 a 

dS"{b) 



d\nb 

dD{b) 



2(5"(6) + 2M(6)(n + 2)/i (3.14) 
= D{b) (3.15) 



din 6 



<'^('')-^^l|F(^+;3 (3.16) 



dlnb \2 OL, 

where n is the number of the field component, /i = /i[T(6), (^"(6), D(6)], /2 = 
/2[T(6), (5" (5), D(6)] and /s = /3[T(6)] are complicated functions but can be approximated in the 
limit of weak dependence of (5" (6) and -D(6) for (5" (6), -D(6) ^ 1. If we are near a QCP which is 
usually at very low temperature the scaling equations will be linearized, keeping only the linear 
term in the coupling constant. The renormalization procedure is stopped at 

(5"(6) = 1 (3.17) 
^From the linearized eqs. (3.12)-(3.14) we get: 

r(6) = Tft^/" (3.18) 
u{b) = «6^/2-2/° (3.19) 
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cen.\ — ^2 In 6 



(5" (6) = e 



(3.20) 



These equations will be analyzed in two cases which are in fact the low temperature and 
high temperature regimes. 

In low temperature limit the integral from eq. (3.20) can be approximated as: 



/= /'"^dxe-="(V2+2/a)j^(yg2r./a) ^ _ 

Jo 1 



/i(0) 



/2 + 2/a 



and eq. (3.20) becomes 



5" (6) = 6^ 
If we introduce the parameter 



^,^2^(n + 2)/i(0) 



1/2 + 2/a 



+ 2Bu{n + 2)ri+"/^ 



2^(n + 2)/i(0) 
1/2 + 2/a 



the low temperature regime is defined by 



(3.21) 



(3.22) 



(3.23) 



T{b) < 1 

and the high temperature regime is defined by 



(3.24) 



T{b) » 1 (3.25) 
^Prom eq. (3.22) using the condition (5" (6) = 1 we get the condition for the low temperature 



regime as 



l/a 

ra 



< 1 



(3.26) 



In the high temperature regime defined now by the inequation (3.26) reversed, we approxi- 
mate the function /i (T) as /i (T) = CT and introduce the new variable v = uT. The linearized 
scaling equations are: 



dv{b) 3 
dl^ = 2^^^) 



(3.27) 



d5"(6) 
din 6 



25'^{b) + 2Cv{h){n + 2) 



(3.28) 
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and the gaussian behavior appears if d°'{b) = 1 and v{b) <^ 1. The general solution of eqs. (3.27) 
and (3.28) have the form 



v{b) = ve 



77^3 In 6/2 



(3.29) 



(5"(6) = e 



2 In; 



6 + 2Cv 



(3.30) 



with C = 2(n + 2)C. 

The initial conditions have been determined following the procedure from Ref. [11] 



T = T- 



ra + 2Bu{n + 2)r^+"/^ 
where B is a constant defined in Ref. [11]. The condition v{h) ^ 1 is defined if: 



(3.31) 



R 



< 1 



(3.32) 



[ra + 2{B + C)u{n + 2)T^+'^/^f'^ 
which is in fact the condition for the validity of the Gaussian model (Ginsburg criterion). The 
critical temperature is well approximated by: 



4/(4+a) 



_2{B + C)u{n + 2)_ 
This result for Tc{5) can be obtained from eq. (3.33) as: 



(3.33) 



T ~ ^4a/(4+a) 



(3.34) 



and for a = 4/5, Tc~ (5° 
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3.3 The Specific Heat 

In order to calculate the specific heat we will use eq. (3.10) for the free energy F. Neglecting in 
the lower approximation the second term we obtain 



F{T) = F(6)6-2/°^-^/2 
The exact solution of eq. (3.16) has the form 



(3.35) 
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= 62/"+5/2 /■'''^^2.g-(2/a+5/2)xyg(2.g2x/a) (ggg) 

Jo 

In order to calculate the temperature dependence of the free energy expressed by eq. (3.36) 
we consider the variable x in the domains 

0<x<|lni (3.37) 

a 1 

-ln-<a;<ln6* (3.38) 

where b* is defined by (5" (6*) = 1 and from eq. (3.17) was calculated as b* = T~"/^. Following 
Ref.[ll] in the first domain /3(r) ^ C3T and in the second domain /3(r) ^ DT. Using these 
approximations we obtain from eq. (3.36) 



^-2/a-5/2 ^(-^^ ^ ^ j.l+5a/4 



C3 dTiT-'"''^ + D dTiT-"""^^ 
where T\ = rexp[2a;/Q;]. If we take a = 4/5 from eq. (3.39) we calculate 



(3.39) 



F{T) = -CaT^ In - + -DT^ - -DTb*~^/^ (3.40) 
5 T 5 5 



and 7(r) = C„/r as: 



7(r)=7eo + 7ln^ + o(^i2) (3-41) 
a result which shows that using RNG for the phenomenological model with a = 4/5 we obtain 
the ln(l/T) term in 7(T)which is in fact done by the non- Fermi behavior of the model. 

Recently Ramazashvilli [13] used the same method studying QLP for such a model with 
a = 1. Our results are consistent with the results from Ref. [13], excepting the specific heat 
coefficient which has T^/^ dependence. The phase diagram can be calculated from the free 
energy 

F ~ ((5" + Dq^ + g^) (3.42) 

and for D > 0, 6°^ we get an ordered phase at g = 0. A second ordered phase can be obtained for 
D < 0, q = ib(£)/2)^ separates the two ordered phases from disordered phase. The importance 
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of the and qj_ terms in the susceptibihty can be also discussed as in Ref. [13] and we get the 
same result, the only difference being that S{b) has to be replace with S°'{b). 

3.4 Summary of the Results 

•We calculated the specific heat dependence of a 3D itinerant electron system near a Lifshitz 
point using the phenomenological model considering the susceptibility of the form ~ + 
S'^ + \lo\°'+T'^. This form describes a similar behavior with the phenomenological model proposed 
in Ref. [7,8] and using the T 7^ RNG proposed by Millis we showed that the specific heat 
presents a InT contribution for a = 4/5 which is typical behavior for a non- Fermi system. 

Our results can be considered as generalization of the recent results obtained in Ref. [13] 
where the exponent a has the value a = 1. The first calculation [11] of the specific heat using 
RNG methods near a QCP with zero critical temperature used for the dynamic exponent z = 
2, ^; = 3 or z = 4 but no InT dependence has been obtained. 
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Chapter 4 



Conclusion 

We conclude this thesis presenting the original results obtained and mentioning the unsolved 
problems in this field. 
Chapter 1: 

(1) The temperature dependence of a marginal Fermi liquid has been calculated . We showed 
that the expected TlnT correction is characteristic for the low temperature domain. The high 
temperature domain has a supplementary correction. The results are in agreement with the 
non- Fermi behavior of some metallic systems in the low temperature domain, (results contained 
in: M.Crisan, CP. Moca Journal of Superconductivity 9 49 (1996)) 

(2) We calculated the self-energy at T = for a two dimensional fermionic system with 
hyperbolic dispersion. The existence of the saddle points in the energy gives rise to a marginal 
behavior, a result which has been obtained by numerical calculations. We present the many- 
body calculation of the sclf-cncrgy. Wc showed that even in a (RPA) approximation this model 
present a marginal Fermi liquid behavior. (results contained in: CP. Moca, M. Crisan Journal 
of superconductivity 10 3 (1997) ) 

(3) In order to give a stronger support to this model we adopted the Renormalization Group 
Method for this model.. The calculation of the wave function renormalization constant Z will be 
performed using the "poor man's renormalization" method and we will show that Z — > with 
the infrared cut-off A as Z{A) = where C is a constant, (results contained in: M. Crisan, C. 
P. Moca Modern Physics Letters B9 1753 (1995)) 

The last experimental results and the theoretical investigation showed that the non-Fermi 
behavior is associated with the two novel features of these systems, namely the occurrence of 
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the pseudogap and the proximity of a Quantum Phase Transition. This problems have 
been studied by different authors for a normal Fermi-liquid. However we adopted the Anderson 
non-Fermi model because it can be reduced to the normal model taking the parameter a equal 
zero. 

Chapter 2 

(4) The electronic self-energy due to the electron-spin interaction is calculated using the 
one-loop approximation for the two dimensional system and quasi-two dimensional (anisotropic) 
model. We analyzed the relevance of the diffusive modes and the temperature dependence of 
the magnetic correlation length for a possible temperature dependence of the pseudogap. The 
results contained in this section are in the spirit of the method presented by the Tremblay group 
for the normal Fermi systems, (results contained in: CP. Moca, I.Tifrea, M. Crisan (accepted 
for publication in Journal of Superconductivity) ) 

(5) We study the influence of the amplitude fluctuations of a non-Fermi superconductor on 
the energy spectrum of the two-dimensional Anderson non-Fermi system. In order to perform 
such a calculation we had to elaborate Schmid self-consistent model of fluctuations in BCS 
superconductor for the fluctuations in an Anderson non-Fermi superconductors. The new prop- 
agator of fluctuations have been calculated and in the limit of a = it gives the results from 
BCS. These fluctuations are also classical and give a temperature dependence in the pseudogap 
induced in the fermionic excitations of the Anderson model, (results contained in: M. Crisan, 
CP. Moca, I. Tifrea Phys.Rev. B59 14680 (1999)) 

In order to study the possibility of the superconducting state in the Anderson model we also 
studied the problem formulated by the Nozieres and Schmitt-Rink who considered the crossover 
between weak (BCS) and strong coupling (BEC) using a Fermi liquid model. 

(6) Using the field-theoretical methods we studied the evolution from BCS description of a 
non- Fermi superconductor to that of Bose-Einstein condensation (BEC) in one loop approxima- 
tion. We showed that the repulsive interaction between composite bosons is determined by the 
exponent a of the Anderson propagator in a two dimensional model. For a ^ the crossover 
is also continuous and for a = we obtain the case of the Fermi liquid, (results contained 
in:I. Tifrea, CP. Moca, M. Crisan (presented at the 6*^ International Conference, Materials 
and Mechanisms of Superconductivity and High Tem,perature Superconductors, February 20-25, 
2000, Houston, Texas, USA , to be published in Physica C) 

The proximity of a Quantum Phase Transition, recently verified by many experimental re- 
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suits, has been studied using the RNG method for the phenomenological model proposed by the 
group of Hilbert von Lohneysen. 
Chapter 3 

(7) Using the renormalization group approach proposed by Millis for the itinerant electron 
systems, in the case oi d = z = 2.5 we calculated the specific heat coefficient 7(T) for the mag- 
netic fluctuations with susceptibility ~ + + fil) near a Lifshitz point. The constant 
value for a = 4/5 and the logarithmic temperature dependence, specific heat for the non- Fermi 
behavior, have been obtained in agreement with the experimental data. (results contained in: 
CP. Moca, I. Tifrea, M. Crisan (accepted for publication in Phys. Rev. B) ) 

Open problems: 

In spite of fact that the non-Fermi model proposed by Varma et.al. and Anderson may 
explain many experimental results the new very accurate measurements showed that there are 
still open problems which have to be explained even in a scmi-phenomenological model. In this 
respect a non-Fermi liquid model have to explain quantitatively the 

• peak-dip-hump structure of the ARPES lineshape and tunnel spectroscopy. 

• the pseudogap spectrum in vortex core of underdoped cuprates superconductors 

• the existence of resonances in neutron scattering 

• marginal or non-Fermi liquid like relaxation rates in optical conductivity 

• nature of the magnetic fluctuations and the magnetic properties of the non-Fermi liquid 
as well as the transport in (pseudo) two-dimensional non-Fermi systems. In this problem we got 

some new results studing the electron-hole channel for the non-Fermi Anderson Model (results 
contained in: M. Crisan, CP. Moca, I. Tifrea (to be published ) ) 
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